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It is widely assumed that a star and its protoplanetary disk are ini- 
tially aligned, with the stellar equator parallel to the disk plane. When 
observations reveal a misalignment between stellar rotation and the 
orbital motion of a planet, the usual interpretation is that the initial 
alignment was upset by gravitational perturbations that took place 
after planet formation. Most of the previously known misalignments 
involve isolated hot Jupiters, for which planet-planet scattering or 
secular effects from a wider-orbiting planet are the leading explana- 
tions. In theory, star/disk misalignments can result from turbulence 
during star formation or the gravitational torque of a wide-orbiting 
companion star, but no definite examples of this scenario are known. 
An ideal example would combine a coplanar system of multiple plan- 
ets — ruling out planet-planet scattering or other disruptive post- 
formation events — with a backward-rotating star, a condition that 
is easier to obtain from a primordial misalignment than from post- 
formation perturbations. There are two previously known examples 
of a misaligned star in a coplanar multi-planet system, but in nei- 
ther case has a suitable companion star been identified, nor is the 
stellar rotation known to be retrograde. Here, we show that the star 
K2-290 A is tilted by 124 + 6 degrees compared to the orbits of both 
of its known planets, and has a wide-orbiting stellar companion that 
is capable of having tilted the protoplanetary disk. The system pro- 
vides the clearest demonstration that stars and protoplanetary disks 
can become grossly misaligned due to the gravitational torque from 
a neighbouring star. 
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formation 


he K2-290 system (1) consists of three stars. The primary 

star, K2-290 A, is a late-type F-star with a mass of 1.19 + 
0.07 solar masses (Mo) and a radius of 1.51 + 0.07 solar radii 
(Ro). The secondary star, K2-290 B, is an M-dwarf with a 
projected orbital separation of 113 + 2 astronomical units (au). 
The tertiary star, K2-290 C, is another M-dwarf located further 
away, with a projected separation of 24671177 au. The primary 
star harbors two transiting planets. The inner planet “b” has 
an orbital period of 9.2 days and a radius of 3.06 + 0.16 Earth 
radii (Rẹ), making it a “hot sub-Neptune.” The outer planet 


c” is a “warm Jupiter” with orbital period 48.4 days, radius 
11.3 +40.6 Re, and mass 246 + 15 Earth masses (Mẹ). 


Stellar obliquity in the K2-290 A system 


We performed high-resolution optical spectroscopy of K2-290 A 
over time intervals spanning the transits of both planets, as a 
way of measuring the sky-projected obliquity of the star via 
the Rossiter-McLaughlin (RM) effect. The physical basis of 
this method is that the transiting planet blocks a portion of 


www.pnas.org/cgi/doi/10.1073/pnas.2017418118 


the stellar photosphere, leading to a distortion in the star’s 
absorption lines. When the planet is in front of the approaching 
half of the rotating star, it blocks some of the light that would 
ordinarily be blueshifted due to the Doppler effect associated 
with stellar rotation. The absorption lines show a small deficit 
on the blue side, which can be detected directly or as an 
apparent redshift of the entire line. The time evolution of 
the spectral distortions throughout a transit depends on the 
sky-projected stellar rotational velocity (vsini,) and the sky 
projection of the stellar obliquity (A). For example, if the 
stellar rotation and the planet’s orbital motion are aligned, 
then the planet blocks the blueshifted (approaching) side of 
the star during the first half of the transit, and the redshifted 
(receding) side of the star during the second half of the transit. 
This would be observed as an anomalous redshift followed by 
a blueshift. 


Planet c. This is the opposite of the pattern that was observed 
for the outer planet of K2-290 (Fig. 1 left panel). The observed 
pattern of radial velocities implies that the planet’s orbit is 
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Fig. 1. The Rossiter-McLaughlin effect for both of the planets in the K2-290 system. Left: Data for the outer, larger planet c. Gray circles are 
radial-velocity (RV) data obtained with HARPS-N on 25 April 2019. White circles show HDS RVs obtained on 12 June 2019. The error bars indicate the internal uncertainties 
as derived for the RV data points by the Data Reduction Software of the spectrographs. The solid gray line is the best-fitting model, which has a projected obliquity of 
Ac = 153 + 8 deg. The lower panel shows Observed minus Calculated (O-C) between the data and the best-fitting model. The dashed gray line is a model in which A« is zero 
and all the other parameters are the same as in the best-fitting model. The horizontal bars denote the time intervals from first to fourth contact (light gray) and second to third 
contact (dark gray). Right: Same as left, but for the smaller, inner planet b. These RVs were obtained on 20 July 2019 with ESPRESSO. We find A, = 173733 deg. The 
shaded gray area shows the time interval when our view of K2-290 was blocked by clouds. 


retrograde with respect to the star’s rotation. The data are 
based on observations of two different transits. We observed 
the first transit on 25 April 2019 with the High Accuracy 
Radial-velocity Planet Searcher North, [HARPS-N, (2)] on the 
3.6m Telescopio Nazionale Galileo on La Palma, in the Canary 
Islands. For the second transit, on 12 June 2019, we used the 
High Dispersion Spectrograph, [HDS, (3)] on the 8.2m Subaru 
Telescope on Mauna Kea, in Hawaii. Together, these two data 
sets provide complete phase coverage of the 8.1-hour transit 
(Supplementary Information). 

By fitting a parameterized model to the RV time series 
(see Supplementary Information), we find the sky-projected 
obliquity to be A. = 153 + 8 deg. We also find vsini, = 
6.9702 km s7! (Fig. 1), which is consistent with the value 
that was determined independently from the observed line- 
broadening of the stellar spectrum (1). 


Planet b. What about the orbital orientation of the inner 
planet? The observation of the RM effect for the inner planet 
is more challenging because of the planet’s smaller size. One 
might expect the orbits of the two planets to be closely aligned, 
based on prior observations and statistical analyses of the Ke- 
pler multi-planet systems (4-6), but K2-290 is unusual in 
having a giant planet. Furthermore, the statistical studies 
could not tell whether the planets always orbit in the same di- 
rection, or if they could sometimes orbit in opposite directions 
(an admittedly speculative possibility). 


Limits on mutual orbital inclinations from long-term dynamical simu- 
lations We performed numerical integrations of the gravita- 
tional dynamics of the two-planet system in order to check on 
the long-term stability of the system for different choices of 
the mutual orbital inclination. We used the mercury6 code 
(7), including the effects of general relativity. We assumed 
both orbits to be initially circular. This is a conservative 
assumption in that sense that any initial eccentricity would 
extend the range of unstable mutual inclinations. The copla- 
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nar configuration is stable for at least 2.8 Gyr, the maximum 
timespan that was simulated. An anti-aligned configuration, 
with a mutual inclination of 180 deg, is also stable over the 
same timespan. 


We tried a configuration in which planet c’s orbit is mis- 
aligned relative to the star by the observed amount, and planet 
b is aligned with the stellar equator. This configuration is also 
stable for at least 2.8 Gyr, although both planets do not often 
transit at the same time. When one planet transits, the other 
planet also transits about 15% of the time. 


Based on these integrations, we calculated the expected 
level of transit-timing variations (TTV) and transit duration 
variations (TDVs) using the code described in Ref. (8). In 
both scenarios, the expected level of TTVs is too small to be 
detectable using the currently available data set. The TT Vs 
are on the order of 0.1 minutes for the inner planet, and 0.01 to 
0.1 minutes for the outer planet (depending on the mass of the 
inner planet). For the case in which the inner planet is aligned 
with the stellar spin, the TDVs of the inner planet are on the 
order of 10 minutes over three years. Our simulations also 
showed that values of the mutual orbital inclinations between 
74° and 112° are unstable, as Lidov-Kozai cycles of planet b’s 
orbit drive up its eccentricity. This makes the planet prone to 
collisions with the star, or tidal disruption. We display these 
limits in Fig. 2. 

We also tested how different values of the planetary masses 
and orbital eccentricities influence these stability zones. The 
most significant influence is from the orbital eccentricity of 
planet c. For example, taking the eccentricity to be 0.144, 
we found that the unstable range of mutual inclinations is 
enlarged to 59° — 136°. The best available constraint on the 
orbital eccentricity based on RV data is a 3-o upper limit of 
0.24 (1). 

In summary, we found that nearly perpendicular configura- 
tions can be ruled out, but mutual inclinations near 0° and 
180° are both viable. 
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Fig. 2. Results for the projected obliquities of K2-290 with re- 
spect to both of the known planets. The upper panel shows the posterior 
for the projected obliquity of planet c, based on the RM measurements shown in 
Fig. 1. The lower panel shows the same for the smaller planet b (black). The orange 
posterior also takes into account the constraints form our orbital stability calculations. 


Spectroscopic transit observations - planet b To decide be- 
tween these possibilities, we observed a transit of the inner 
planet with the newly commissioned Echelle Spectrograph 
for Rocky Exoplanets and Stable Spectroscopic Observations 
[ESPRESSO, (9)] and one of the 8.4m Very Large Telescopes 
at Paranal Observatory, in Chile (Supplementary Information). 
The data, shown in the right panel of Fig. 1, right panel, rule 
out anti-aligned orbits and are consistent with aligned orbits 
(Fig. 2). For the determination of the projected obliquity 
of planet b, we followed the same approach as for planet c, 
this time setting priors on vsin i, and the Kepler-band limb- 
darkening parameters based on the previous analysis. We used 
the same constraints on the limb darkening coefficients for the 
ESPRESSO data that were used for the HARPS-N data. 

The result is Ap = 173° 5. deg (Supplementary Materials), 
which is compatible with the value obtained for the outer 
planet. The lower signal-to-noise ratio means we cannot mea- 
sure Àp as well as Ac. Indeed, if we did not already have 
precise prior knowledge of the transit depth, the transit im- 
pact parameter, the mid-transit time, and vsin ią, we could 
not have been confident that the RM effect for planet b was 
detected. Our prior knowledge of those parameters guarantees 
that the RV anomaly has an amplitude of a few meters per 
second (10). The only effectively free parameters in the model 
are Ap and the velocity zeropoint. 

For an even simpler test, we fitted the data with two 
different models: one in which the two planets are exactly 
aligned (Ab = àc), and the other in which they are anti-aligned 
(Ap = 180 deg + àc). There are no free parameters in either 
model apart from the velocity zeropoint. The aligned model 
has a x” statistic of 13, while the anti-aligned model is a 
poorer fit with x? = 26. In both cases, the number of degrees 
of freedom is 13. Based on these tests, we conclude that the 
data favor the aligned case, although they do allow a small 
probability (on the order of a few percent) for the even more 
exotic possibility of anti-aligned orbits. 

As a consistency check, we also investigated if the RM effect 
can be detected as deformations of the stellar lines during 
the transit (i.e. the Doppler shadow — see Supplementary 
Information). We found: (i) The “planet shadow” cannot be 
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detected by eye in the time series of cross-correlation functions. 
(ii) This visual non-detection is consistent with the expected 
strength of the deformation and the Signal-to-Noise Ratio 
(SNR) of the data (Fig. S4, panels A-C). (iii) By stacking 
all the ESPRESSO observations we detected a signal with 
the expected characteristics, from which we measured Ay = 
157+34°, consistent with the RV-based result and the expected 
strength of the RM/planet shadow signal (Fig. S4, panels D-F). 
We have more confidence in the analysis of the anomalous 
RVs, because the method is simpler than the analysis of the 
Doppler shadow. 


Stellar rotation period and stellar inclination. We determined 
the stellar rotation period based on the quasi-periodic bright- 
ness fluctuations that are seen in the 80-day time series ob- 
served with the Kepler telescope. For this determination, we 
began with the “K2SFF” version of the photometric time series 
for K2-290 A (11). We removed an overall downward-sloping 
trend by applying a high-pass median filter with a width of 
1,000 time samples (20.8 days). To prepare to compute the 
autocorrelation function, we created a light curve with uniform 
time sampling by filling in the gaps with linear interpolation. 
The resulting light curve is shown in the top panel of Fig. 3. 
The middle panel shows the autocorrelation strength as a func- 
tion of lag, which displays a series of regularly spaced peaks. 
By fitting a linear function to the measured peak location 
versus cycle number (bottom panel), we found the slope to 
be 6.63 + 0.06 days. We adopted this value for the rotation 
period of K2-290 A, although we enlarged the fractional uncer- 
tainty from 1% to 10% to account for systematic errors due 
to differential rotation (12, 13). 

Based upon this value, we calculated the equatorial rotation 
velocity v under the assumption of uniform rotation: 


=11.2+ 1.1 km s™*, (1] 


where we have used R, = 1.511 40.075 Ro (1). This value 
for the rotation velocity is larger than the measured projected 
rotation velocity, vsini, = 6.9+0.5 km s~', implying that 
sini, is less than unity. Since sin iorb is known to be close to 
unity because of the detection of transits, we have here an 
independent line of evidence for a large stellar obliquity. 

We determined the inclination angle of the stellar rotation 
axis, and the obliquity of the star, Y, using a variation of the 
Bayesian inference technique advocated by (14). The model 
parameters were R,, Prot, cos ix, for which uniform priors were 
adopted. The likelihood function was taken to be 


2 
ee (a) 


0.075 


Prot —6.63d\? | (vu—6.9km/s\? 
i 1d ) l 0.5 km/s » PI 
where v = 27 R,/Prot and u = y1 — cos? i4. The result was 
sini, = 0.63 + 0.09. Thus we found independent evidence for 
a large stellar obliquity based on the combination of the stellar 
rotation period, radius, and projected rotation velocity. There 
are two possible solutions for 74, which are near 39° and 141°. 
Combining the posteriors of these two solutions with equal 
weight, and including the posterior for A» based on the RM 
effect, the result for the stellar obliquity is 124 + 6 degrees 
(Table 1). 
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Fig. 3. Determination of the stellar rotation period. Top: Detrended 
Kepler light curve, including linear interpolation across the gaps in the time series. 
Middle: Autocorrelation as a function of lag. Bottom: Locations of the sequence of 
autocorrelation peaks in the autocorrelation function (ACF), along with a linear fit. The 
slope of this line is taken to be the rotation period of K2-290A. 


Radial velocity limits on wide-orbiting bodies. By combining 
two new radial-velocity observations taken with HARPS-N 
in summer 2019 (the last two entries of Tab. S1) with the 
HARPS-N measurements obtained previously (1), we have 
data extending over 500 days. We modeled the planetary 
signals along with a possible long-term constant acceleration, 
finding the acceleration to be y =9+5 ms‘ yr™t. Based on 
the estimated mass and orbital distance of K2-290 B (0.368 4 
0.021 Mo and 113 + 2 au), we would expect to see a radial 
acceleration on the order of 5 ms~' yr~'. Therefore, the 
observed long-term acceleration is compatible with zero (within 
2-c) and is also compatible with the expected contribution 
from the nearby M-dwarf. There is no evidence for any other 
wide-orbiting bodies. 


Discussion 


Systems similar to K2-290, with coplanar planetary orbits and 
a grossly misaligned host star (Fig. 4), had been predicted 
to exist as a consequence of the tidal torque on a protoplane- 
tary disk from a neighboring star (15, 16). Another possible 
explanation for such systems is the tilting torque exerted on 
the inner system of planets by a massive planet on a wide 
and highly inclined orbit. The Kepler-56 system features two 
planets on coplanar orbits and a misaligned star (17), and in 
that case, a wider-orbiting third planet was detected through 
long-term radial-velocity monitoring (18). Based on the mass 
and orbital distance of the third planet, it is possible or even 
probable that the planet was responsible for tilting the orbital 
plane of the inner two planets long after the planets formed 
(19). Likewise, the HD 3167 multi-planet system was recently 
found to have a misaligned star (20), but there is not yet any 
evidence for either a wider-orbiting planet or a companion 
star. Turbulence (21) and disk-torquing (22) can lead to mis- 
aligned protoplanetary disks. However retrograde orbits, as 
observed for K2-290A, are difficult to achieve via turbulence 
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Table 1. Selected parameters of the K2-290 system. 


Parameter Value 
Stellar rotation period, P,ot (days) 6.63 + 0.66 
Stellar inclination angle, i, (deg) 3947 
Projected stellar rotation velocity, v sini, (km s71) 6.90.5 
Projected obliquity with respect to planet b, A» (deg) 173155 
Projected obliquity with respect to planet c, A. (deg) 153+ 8 
Obliquity with respect to planet c, Ye (deg) 124+6 


and late infall of material will lead to a further reduction of 
any misalignment (23). 

The unique aspect of K2-290 is that a companion star 
has been detected (K2-290 B) with properties that make it a 
good candidate for the misalignment of the protoplanetary 
disk. There is no evidence for a wider-orbiting massive planet: 
the upper limit on any long-term radial acceleration is about 
10 times smaller than the acceleration that was observed for 
Kepler-56. In addition, star/disk misalignment is an attrac- 
tive explanation for K2-290 because it can easily produce 
retrograde orbits (24, 25). This is because the orientation 
of the orbital plane of a wide binary star may have only a 
weak correlation with the orientation of the protoplanetary 
disk around either star (26, 27). In contrast, producing a 
retrograde system through the action of a wide-orbiting planet 
requires invoking an unseen third planet with an unusually 
high orbital inclination which would be difficult to achieve 
through planet-planet interactions (28). 


Primordial disk/star misalignment scenario. The originally 
proposed mechanism for star/disk misalignment was nodal 
precession of the disk around the angular momentum vector 
of the binary orbit (15). It was later recognized that the grav- 
itational coupling between the star and disk is also important, 
and that misalignments are more likely to arise from secular 
resonances between spin and nodal precession (16). Further- 
more, the magnetic fields of stars as massive as K2-290 A are 
probably too weak to have enforced star-disk alignment (29). 

To demonstrate that this scenario is plausible for the case 
of K2-290, we calculated the system’s secular evolution due to 
mutual gravitational torques between the rotational bulge of 
the host star, the protoplanetary disk, and the companion star 
K2-290 B. We employed the model described in Ref. (25). In 
this model, there is a star of mass M,, radius R,, and rotation 
frequency 2,. To represent the pre-main sequence phase of 
stellar evolution, the stellar radius is set equal to R, = 2Roe 
and the rotation rate is set such that Q, = 0,/.,/GM,/R3 = 
0.1. Including contraction of the stellar radius, and evolution 
of the stellar spin during the disk-hosting phase, does not have 
a significant impact on the ensuing star-disk-binary dynamics 
(16). The star is surrounded by a circular flat disk with an 
inner radius of rin = 4 R,, and an outer radius of rout = 50 au. 
The disk’s surface density profile is 


Ma (t) 


(r,t) =~ ——— 
(r, ) 2 Toutr” [3] 
where 
Malt) = eee 4 


is the disk mass, ty = 0.5 Myr is the viscous timecale, and 
Mao = 0.1 Mo is the initial disk mass. Although modifying the 
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disk properties does modify the likelihood of secular resonance 
crossing, we find a secular resonance occurs over a wide swath 
of reasonable disk parameters for this system (Supplementary 
Materials). 

The model assumes that the two planets form within the 
disk at the locations we observe them today. Our model does 
not take into account the effects of planet migration or photo- 
ionization of the disk as previous work has shown that these 
effects tend to lead to even greater excitation of the stellar 
obliquity during the disk-hosting phase (25). The mass of the 
outer planet is set equal to its currently observed value, while 
the mass of planet b is set equal to 21.1 Mẹ, the 3-0 upper 
limit that was derived from radial-velocity observations (1). 
In general, smaller values of M» increase the chance of large 
star-disk misalignments occurring. The model also includes a 
companion star of mass Mp in a circular and inclined orbit 
with radius ap. We assumed the binary’s semi-major axis is 
greater than the observed projected separation [ap > 113au 
(1)]. We used the secular equations from (25) for the dynamical 
evolution of planet-forming star-disk-binary systems. To these, 
we added the gravitational influence of planets b and c on the 
star-disk ©, and disk-star ©), precession frequencies [see eqs. 
65-66 of (25)]: 


wa = Wsa + Wsb + Wsc, [5] 


Õhs = Was + (Lr /La)wbs + (Le/La)wes, [6] 


where La œ (2/3)MaVGM, rout is the total disk orbital an- 
gular momentum, Li = MiV GMa; are the planets’ orbital 
angular momenta (where i is either b or c), and the precession 
frequencies are 


_ 3kg (2) y GM, R3 (7 


Wsi * 
Oke M, a3 i 


2 
a= Sega (2) 
2 Ai 


[8 


In Eq. (5) and Eq. (6), we use the notation of (25), where 
precession frequencies with tildes are averaged (integrated 
over the radial extent of the disk. We assumed the primary 
star to have moment-of-inertia constants of k, = 0.2 overall 
and ką = 0.1 for the rotational bulge, as appropriate for the 
pre-main sequence phase (24). We neglected the torque on 
the planets from star B (@4p ~ @ap). 

Fig. 5 displays the main result: the time evolution and 
excitation of primordial misalignments or mutual star-disk 
inclinations sa = cos™t (8 . ia) excited by the binary com- 
panion, for different initial disk-binary mutual inclinations 
Ba» = cos*{1a(0) - fg] (where ŝ is the host star’s stellar spin 
axis, ia is the disk’s orbital angular momentum unit vector, 
ip is the binary’s orbital angular momentum unit vector) and 
binary semi-major axis values ag. After a secular resonance 
occurs in the system (@1, ~ @ap), large primordial misalign- 
ments are generated, which are consistent with the measured 
stellar obliquity Ye = 124+ 6 deg over a range of parameter 
values. Since we took the inner planet’s mass to be the 3-0 
upper limit obtained from RV measurements (1), the sup- 
pression of primordial misalignments by short-period, massive 
planets would not occur for this system (25). A misalignment 
would only be averted if the disk dissipation timescale were 
extremely fast (< 0.05 Myr) or if the disk were extremely 
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Fig. 4. Illustration of the architecture of K2-290 A and its two 
known planets. The view is from above the system’s ecliptic plane, along the 
orbital north pole of the two planets. The ratio of the sizes of the orbits is accurate. 
Relative to the scale factor of the orbits, the star has been enlarged by a factor of 5, 
and the planets have been enlarged by a factor of 500. The blue contours indicate 
the 1 and 2 — o confidence intervals for the stellar obliquity, with respect to planet c. 
The red point represents the stellar south pole. The red arrows indicate the sense of 
stellar rotation, and the blue arrows indicate the sense of orbital motion. Our actual 
view of this system is from the side, allowing for the observation of transits. 


compact [< 3 au, (25)], and such disks have to our knowledge 
not been observed. 


For a wide range of possibilities for the angle between 
the disk and the binary orbit, and for a plausible range of 
values for the binary orbital distance, the final value of the 
stellar obliquity can be made consistent with the measurements 
(Fig. 5). The main requirement for producing a large star/disk 
misalignment is the occurrence of a secular resonance. In our 
model, a one-to-one resonance occurs between the precession 
frequency of the stellar spin axis around the disk’s angular 
momentum vector, and the precession frequency of the disk’s 
angular momentum vector around the angular momentum of 
the binary orbit. Previous work found that secular resonance 
crossing occurs as long as the disk loses its mass over timescales 
between 0.05 and 1 Myr (16, 25), consistent with measurements 
of protostellar disk masses and host star accretion rates in the 
Lupus cluster (30). 


In a subset of the models we considered, the companion 
star continues to reorient the orbital plane of the two planets 
after the gaseous disk disappears (Supplementary Materials), 
through nodal precession or another resonance crossing, but 
in all such cases the original misalignment occurs while the 
gaseous disk is still present (Fig. S7). We also note that while 
K2-290 has retained a stellar companion, many systems lose 
their companions after a few Myr (31). Therefore the disk- 
torquing mechanism described here might have operated even 
in systems which are now observed as single stars. 


In short, we find the primordial disk/star misalignment 
scenario to provide a plausible explanation for the observed 
properties of K2-290. It will be interesting to perform similar 
measurements of other systems to see how frequently this 
scenario occurs, keeping in mind that there may be more than 
one way to misalign a disk. Studies of directly imaged disks 
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Fig. 5. Simulation of primordial misalignment from star-disk- 
binary interactions (with embedded planets). The top panels show 
the time evolution of the angle between the angular momentum vectors of the host 
star and the disk (0,4), for different values of the initial disk-binary inclination 045 (0) 
and binary semi-major axis (ap), as indicated in the legend. In all 3 cases shown, 
sa periodically crosses the interval corresponding to the observed stellar obliquity 
(dashed black line with gray band representing the uncertainty). The bottom panel 
shows the evolution of the star-disk @/, and disk-binary ap precession frequencies. 
Large primordial misalignments 6,4 are excited when secular resonance occurs in 
the protostellar system (4, ~ @ap). 


have not yet turned up any examples of a misaligned star 
however, there are hints, albeit limited, from (32) who found 
evidence that 1/3 of their sample had non-aligned configura- 
tions. There are some cases in which the inner and outer disk 
are misaligned, and obviously the star must be misaligned 
with one of them (33). For now, the architecture of the K2-290 
system shows that we cannot safely assume that stars and 
their protoplanetary disks are always well aligned. This may 
lead to a re-interpretation of other systems with observed 
misalignments, especially the hot Jupiters (34). It also helps 
to set expectations for future explorations of exoplanetary 
systems with wider-orbiting planets, more similar to the Solar 
System. 


Materials and Methods 


In this article we analyze high resolution spectroscopic data 
obtained for the K2-290 system during transits of two planets in 
front of their host star. These observations are detailed in the SI, 
section 1. We use the apparent RVs, derived from these observations, 
to determine the projected obliquity of K2-290. We detail our 
choice of model, parameters and prior information in section 2 
of the SI. Here we also describe tests to confirm that our results 
do not critically depend on the exact RV model setup and prior 
information supplied. To supplement these results further we have 
employed alternative analysis techniques for the spectroscopic data 
(sections 3 & 4 of the SI). These techniques do not make use of 
the anomalous RVs observed during transit, rather they rely on 
information contained in the deformation of the stellar line shapes 
during planetary transits. Our analysis of these deformations is 
consistent with our earlier analysis of the RVs. Section 5, the last 
section of the SI, details the setup of our calculations which explore 
possible post-formation secular resonance scenarios. 


The HARPS-N data reported in this paper are archived in the 
INAF Science Archive (https://www.ia2.inaf.it/) under program 
ID A39TAC_ 2. The HDS data reported in this paper are archived in 
the SMOKA Archive (https://smoka.nao.ac.jp) under program 
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ID $19A122. The ESPRESSO data reported in this paper are 
archived in the ESO Science Archive (http://archive.eso.org) 
under program ID 2103.C-5041(A). The Kepler light curve of 
the K2-290 system from the K2 mission reported in this pa- 
per is archived in Mikulski Archive for Space Telescopes, MAST 
(https://archive.stsci.edu/). 

This paper also includes data collected by the K2 mission, 
which was funded by the NASA Science Mission directorate. 
This research made use of Lightkurve, a Python package for 
Kepler and TESS data analysis (Lightkurve Collaboration, 
2018). 
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Supporting Information Text 
1. Spectroscopic transit data 


A. Spectroscopic transit observations of planet c. We observed the first transit on 25 April 2019 with the High Accuracy 
Radial-velocity Planet Searcher North, [HARPS-N, (1)] on the 3.6m Telescopio Nazionale Galileo on La Palma, in the Canary 
Islands. The observations began just before mid-transit and lasted throughout the second half of the transit. There was 
a two-hour interruption when the telescope was commandeered by other observers pursuing a Target of Opportunity. The 
HARPS-N spectra cover the wavelength range from 383 to 693 nm, with a spectral resolution of 115,000. The exposure time 
ranged from 20 to 45 minutes. All HARPS-N observations obtained during the transit night are presented in Tab. S1. 

For the second transit, on 12 June 2019, we used the High Dispersion Spectrograph, [HDS, (2)] on the 8.2m Subaru Telescope 
on Mauna Kea, in Hawaii. With HDS, we obtained 3 reference spectra of K2-290 A on 11 June 2019, the night before a transit 
of the outer planet. These served to establish a radial-velocity baseline. On the transit night of 12 June 2019, we began 
observing at ingress and were able to continue for about 3/4 of the transit duration. We also obtained three pre-ingress spectra, 
but none of them proved useful: the first two spectra were excessively contaminated by twilight, and the third exposure suffered 
from a problem with the detector readout process. For all of the HDS observations, we used Image Slicer #2 and the Iodine 
(I2) absorption cell in the light path. The HDS spectra cover the wavelength range from 494 to 759 nm, although the iodine 
absorption lines are almost all within the range from about 500 to 600 nm. The spectral resolution was approximately 80,000 
and the exposure time ranged from 15 to 20 minutes. Rather than investing Subaru observing time to obtain an iodine-free 
template spectrum, we created a template by co-adding archival HARPS and HARPS-N spectra (which have a higher spectral 
resolution than our HDS spectra). The HDS RVs are listed in Tab. S2. Together, these two data sets provide complete phase 
coverage of the 8.1-hour transit. 


B. Spectroscopic transit observations of planet b. The observation of a transit of the inner planet was performed on 20 July 
2019 with ESPRESSO (3), which was connected to Unit Telescope 3. The ESPRESSO spectra cover the wavelength range 
from 380 to 788 nm, with a spectral resolution of 140,000. The exposure time ranged from 15 to 30 minutes. We extracted 
Cross Correlation Functions (CCFs) and RVs using the standard Data Reduction Software. 

Observations began prior to the ingress and lasted for most of the transit. Clouds began covering the sky at around the 
time of the transit midpoint. The effects were to lower the signal-to-noise ratio of each spectrum, increase the contamination of 
the spectra by moonlight, and increase the strength of telluric spectral features. We chose not to include the affected spectra in 
our analysis, although for completeness, Tab. S3 provides the RVs that were extracted from the affected spectra (marked with 
a x symbol). 


2. Details of determining projected obliquities via anomalous transit RVs 


To measure projected obliquities of K2-290 A with respect to planets b and c, we fitted a parameterized model to time series of 
apparent radial velocities and the photometric time series obtained with the NASA Kepler telescope during the K2 mission. 
The model was based on the one that was described in Ref. (4), extended to include the Rossiter-McLaughlin effect (5). The 
photometric time series was obtained via the Lightkurve tool (6), which we used to extract a light curve from the K2 data. 
The resulting light curve has slightly better noise statistics than the time series used in the discovery paper. The RV time 
series include observations from the transit nights, as well as the three HDS data points obtained the night before a transit of 
planet c. 

The model parameters were the projected spin-orbit angles of the two planets (Ap and Ac), the projected stellar rotation 
velocity (vsini,), quadratic limb darkening parameters in the V-band (u1,v and u2z,v), and terms representing the stellar 
microturbulent and macroturbulent velocities (Umic and Umac). There were also some nuisance parameters relating to the HDS 
and ESPRESSO spectrographs: the additive radial-velocity offsets (“zeropoints”) to align them with the HARPS-N velocity 
scale (yups and yrsp). 

The parameters governed mainly by the photometry were the orbital periods (P,,-), the times of midtransit (Tp,-), the scaled 
orbital separations (ap,-/R,), the radius ratios (rp,-/R.), the orbital inclinations (i,b,c) and the quadratic limb darkening 
parameters (u1,k and u2,K) for the Kepler bandpass. 

We used the best estimate of the stellar density (p, = 0.48510374 gcm™?) from (4) as a prior constraint, which sharpened 
the determination a,,-/R, and the orbital inclinations. The results for the projected obliquities do not depend critically on 
this constraint. Additional prior constraints were placed on the orbital RV semi-amplitude for planet c (Ke = 38.4 + 1.7 ms~°) 
and the 3 — o upper limit of K, < 6.6 ms~* on planet b (4). 

For determining the projected obliquity of planet c we imposed a Gaussian prior constraint on vsini,, with a mean of 
6.5 kms“! and a standard deviation of 1.0 km s~', based on an analysis of the spectral line widths (4). For the limb 
darkening parameters, we used Gaussian priors based on the tabulated values u1,x = 0.38 £0.10 and u2x = 0.14+ 0.10, 
u1 HARPS-N = 0.37 + 0.10 and u2, HARPS-N = 0.31 + 0.10, and ui1,Hps = 0.41 + 0.10 and u2,Hps = 0.33 + 0.10 (7). For the rest 
of the parameters, we used uniform priors. 

To calculate the posterior distribution of the parameters, we used the Markov Chain Monte Carlo Method (MCMC) as 
implemented in the emcee code (8). We used 100 walkers in both runs and let the code running until the chains had reached a 
length of 50 times the autocorrelation time at which point we considered the runs to be converged. For planet c this took 
12,200 draws and 56,300 for planet b. We discarded 2 times the autocorrelation time of the first draws in both runs. 
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In Tab. S4 we report the results for all parameters based on the 15.9%, 50%, and 84.1% levels of the cumulative marginalized 
posterior distributions. Fig. 1 in the main article and Fig. S1 show the RV and phase folded photometry data as well as the 
maximum likelihood model. The two panels in Fig. S2 display the correlation between the projected obliquities and vsin ix. 


We also tested the robustness of A. against the orbital eccentricity of planet c. The authors of the discovery paper (4) prefer 
a circular solution for the orbit, based on the available RVs and photometry. For such a circular orbit we obtain Ae = 153+8 deg. 
If we use their eccentric orbital solution (e = 0.144 + 0.033 and w = 70 +9 deg), then we find A = 156*8 deg. Our measurement 
of the projected obliquity does not strongly depend on the assumed orbital shape. We follow here the discovery paper and 
adopt the circular solution. 


The K2 mission only observed two transits of planet c and while we use information from out of transit RVs covering a time 
period of more than 500 days, RVs are not as powerful in constraining transit timing as transit photometry. To what extend 
does the uncertainty in To and P. and consequently the uncertainty in the particular mid-transit times of our spectroscopic 
observations with HARPS-N and HDS affect Ac? From the posterior of our MCMC run we find the uncertainty interval for the 
particular mid transit timing of our HARPS-N observations to be enlarged to ~10 min, TOHARPs—N = 2458599.578 + 0.005 BJD. 
This can be compared to the To interval listed in Tab. S4 (2458019.1723 + 0.0003 BJD), which reports the timing of twelve 
epochs earlier and closer to the K2 observations. The HDS transit observations had been carried out exactly one epoch after 
the HARPS-N observations and the timing uncertainty for that particular transit is unchanged to the uncertainty for the 
HARPS-N transit. The corner plot in Fig. S3 illustrates that this timing uncertainty does lead to an enlarged uncertainty 
interval for Ae and is accounted for in our final result for A. = 153 + 8 deg. 

Finally we tested if the two RV data sets from HDS and HARPS-N give consistent results independent from each other. For 
this we performed two additional MCMC runs, both identical to the above mentioned run, but with differing RV data sets. In 
the first run we only included HDS RVs and we find Ac = 14475 deg. In the second run we only included the HARPS-N RVs 
and find À = 167734 deg. The two RV data sets obtained with two different instruments during two different transit events 
give consistent results. 


3. Alternative RM analysis for planet b 


As a further check on our results on the projected obliquity of planet b we analyzed the spectra at the level of spectral-line 
distortions (the planet’s “Doppler shadow”), rather than the level of anomalous radial velocities (9-14). For this exercise, we 
used the CCFs generated by the standard ESPRESSO data reduction software. First, we created a mean out-of-transit CCF 
based on the pre-ingress data. Then, we subtracted the mean out-of-transit CCF from the CCF of each ESPRESSO spectrum. 
Fig. S5, panel C, shows these CCF residuals. Any distortion of the CCFs during the transit is not obvious to the eye. To 
proceed, we used a stacking technique (12, 15). For a given choice of Ap and vsin ią, we identified those portions of the CCF 
residuals that should have been affected by the planetary shadow. These residuals were added — weighting them by a function 
representing the local stellar surface velocity field governed by macro-turbulence — to form a statistic that should be largest 
when Ap and vsin i, are chosen correctly. We further weighted this statistic by a Gaussian prior on vsin?, with a mean of 
6.9 km s~! and a standard deviation of 0.6 km s~+, based on the analysis of the RM effect for planet c. Fig. S5, panel F, shows 
the strength of this statistic as a function of Ap and vsin iẹ. 


To help interpret the results, we created simulated data affected by the Doppler shadow, and then analyzed the simulated 
data with the same procedures used on the real data. An important parameter for the simulation is the macro-turbulent velocity: 
the stronger the macro-turbulence, the more the planetary-induced signal is smeared out in wavelength (16). (Micro-turbulence 
and instrumental broadening also contribute to the smearing, but for an F star such as K2-290 A, macro-turbulence is expected 
to be the largest contributor.) Our combined analysis of the real data gave an estimate for the macro-turbulent velocity of 
4.8 + 1.4 km s™+. This is consistent with the value that was obtained by analysing the RM effect for another star of a similar 
effective temperature (11). The scaling relations given in Refs. (17, 18) lead to estimates of 5.4 km s~' and 7.3 km s~’ for 
K2-290 A. Based on all this information, we chose 4.8 km s7} as the basis for our simulations. We constructed simulated 
CCFs for various choices of A and vsin ią, convolved them with a model for the spectrograph broadening function (Fig. S4, 
panel A), and added Gaussian noise to match the signal-to-noise ratio of the data (Fig. S4, panel B). Through this numerical 
experimentation, we concluded that the apparent absence of a signal in Fig. S4, panel C, is consistent with our simulations of 
the shadow given the system parameters and the achieved SNR in the observations. We also carried out the same stacking 
for our simulations with and without noise (Fig. S4, panels D-E) and found the simulations with noise to resemble the weak 
detection of the RM signal shown in Fig. $4, panel F. 


We conducted a bootstrap analysis, creating simulated datasets by drawing randomly from the different CCFs, with 
repetitions allowed, and sampling form the posteriors of all the relevant parameters (period, mid-transit time, orbital inclination, 
scaled orbital separation, stellar projected rotation speed, radius ratio, and macro-turbulence). For the macro-turbulence we 
drew from a normal distribution with mean 4.8 and standard deviation 2 kms~'. For each draw, we calculated the stacked 
signal for different values of Ay. To this stacking signal as function of projected obliquity, we fitted a Gaussian function. We 
repeated this procedure 10,000 times. From the histogram we obtained >» = 157 + 34°. Thus, the analysis of the Doppler 
shadow and the analysis of the anomalous RVs led to similar results. We place more confidence in the analysis of the anomalous 
RVs, because it is simpler than the analysis of the Doppler shadow. 
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4. Alternative RM analysis for planet c 


Fig. 1 in the main article does demonstrate that the RM effect of planet c is clearly detected in RV space. However the plot 
also shows that no pre-ingress or post-egress data had been obtained during the transit nights. Can we confirm the detection 
of the RM effect using a different diagnostic? The HARPS-N spectra do not have iodine absorption lines imprinted. (The 
HDS spectra have been obtained through a iodine absorption line for wavelength calibration.) We therefore focus on the CCFs 
derived from the HARPS-N spectra by the Data Reduction Software of this instrument. We do not perform a classical analysis 
of the "Doppler shadow". The SNR in the CCFs is relatively low compared to the expected signal (Fig. S5) and the observations 
had been interrupted around mid-transit. For these reasons we took the mean of the first 4 CCFs taken before mid-transit and 
subtract this mean from the average of the last 5 transit observations taken after mid-transit. (For a better estimation of the 
scale of the line deformations we have scaled all CCFs by the same factor so that they scale between [0-1]. This results in 
the same scaling as applied in Fig. S4 and the line deformations as caused by the two planets can be compared.) In Fig. S5 
we show next to the data also the expected signal for our best fitting RV model with A. = 153 deg and an aligned model 
Ac = 0 deg. The observed signal is fully consistent with the retrograde orbit. As for planet b we prefer the analysis of the RVs. 


5. Post-formation secular resonance scenario 


In the system’s current configuration, planet b and c are too tightly coupled to the rotating primary star to experience 
significant misalignment by the M-dwarf. However, the system may have encountered a secular spin-orbit resonance after 
the disappearance of the gaseous disk. We investigated this possibility by numerically integrating the differential equations 
representing the coupled evolution of the orbital angular momentum of each planet, the orbital angular momentum of the 
binary orbit, and the stellar angular momentum. For details of this method, see (19). The resonance requires a sufficiently low 
mass for planet b, no more than about 5 Mg, to avoid coupling the planets too strongly to the stellar spin. It also requires 
that the binary orbital distance is not much larger than about 100 au. We assumed the primary star has a moment-of-inertia 
constant of k, = 0.06 overall, and kg = 0.01 for the rotational bulge (19). We further assumed M» = 5 Mg, and adopted a 
somewhat arbitrary initial mutual inclination of 50 deg for the binary orbit relative to the planetary orbital plane. The planets 
were initially aligned with each other and with the stellar spin. 

In the first simulation (Fig. S6, top panel), the star has not yet begun expanding and has an assumed radius of 1.3 Ro, and 
the binary semi-major axis is ag = 100 au. In the second simulation (Fig. S6, bottom panel), the star has its current radius of 
1.511 Ro and we set ag = 80 au, assuming a moderately eccentric orbit. In both cases, the planets remain strongly coupled and 
coplanar to each other, and develop a major misalignment relative to the stellar spin. The planets’ strong gravitational coupling 
relies on planet c’s large mass and the host star’s moderate rotation period: a less-massive exterior planet (M. ~ few Mg), or 
a faster-rotating host star (P-ot ~ 1 day), would not enforce coplanarity under the influence of a misaligned stellar spin or 
companion (20). Based on Eqn. (16) of (20), the initial rotation period must have exceeded ~ 0.3|i, — i.|~1/? day for the two 
planets to have remained aligned at early times. Moreover, if the masses of the planets were switched (Mp ~ Mjup, Mc ~ Mo), 
secular resonances may destabilize the planetary system (21). This scenario does not require extreme fine tuning because the 
system could cross through the secular resonance when the gas disk disappears or as the stellar radius evolves. However these 
same choices of parameters would also have led to disk misalignment at earlier times (Fig. S7). 
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Fig. S1. K2 photometry of K2-290. 
Left: Phase-folded light curve of the transits of planet c. The gray line is the best-fitting model. The lower panel displays 
Observed minus Calculated (O-C). The single point with an error bar is not a data point; it illustrates the 1-o uncertainty 
assigned to each data point. Right: Same, but for planet b. 
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Fig. S2. Posteriors of vsini,, Ac, and Ap. Left: Posterior for the projected obliquity with respect to planet c plotted against the projected stellar rotation. The gray dot 
denote the maximum likelihood values used to depict the models in Fig. 1 and Fig. S1. Right: Same, but for planet b. Note that the v sin i, posterior from the run for planet c 
was used as a prior for planet b. 
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Fig. S3. Posteriors of ATO and A¿. Corner plot for the particular mid transit timing uncertainty for planet c during the night of our HARPS-N observation and for Ac. 
For this transit event we found a posterior of To, HARPS—N = 2458599.578 + 0.005 BUD. For creating this plot we subtracted from the y-axis the mean value and expressed 


the time axis in minutes. 
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Fig. S4. The Doppler shadow of K2-290 Ab in the ESPRESSO data. Panel A: Noiseless simulation of the distortion in the cross-correlation functions 
based on the parameters of the best-fitting model to the light curve and the ESPRESSO RV data (Table S1). The vertical dashed lines mark v sin ią, while the horizontal solid 
and dotted lines mark the times of first and second contact. B: Same, after adding Gaussian random noise with an amplitude set equal to the observed noise level. C: Real data. 
The red outlines in panels A-C mark the positions of particular intensity levels in the simulation. These levels are also indicated in the color bars. The scales are calibrated so 


that the out of transit cross-correlation function has a height of unity. Panels D-F display the strength of the “stacked signal” for different choices of v sini, and Az in the 
simulation (D), the simulation with noise (E), and the real data (F). 
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Fig. S5. Illustration of the line deformations in the HARPS-N data taken during transit of planet c. The solid line with dots indicates the 
average post mid-transit CCF from which the mean pre mid-transit CCF was subtracted. The black-dark gray dashed line indicates the model perdition based on the best fitting 
values from a fit to the photometry and RVs (A. = 153 deg). The back-light gray dashed line shows a similar model with the difference that a well aligned orbit (A. = 0 deg) is 
assumed, a poor representation of the data. We note that the DRS (Data Reduction Software) of HARPS-N delivers CCF datapoints which are separated by 0.25 km s7! in 
velocity space. Given the spectral resolution of HARPS-N, R œ 110 000, the FWHM (Full Width Half Maximum) of the PSF (Point Spread Function) is + 2.72 kms~+. This 
leads to an oversampling of the CCF and therefore a "smooth" appearance of the noise. Adjacent datapoints are not completely independent from each other. 
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Fig. S6. Misalignment generated by secular resonance. Mutual inclination between the host star’s equatorial plane and planet c’s orbital plane (sc) for two 
combinations of M-dwarf separation and primary star radius. 
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Fig. S7. Primordial and post-formation misalignment parameter space. Stellar rotation periods P»ot and re-scaled semi-major axis @g = ag (1 — 
e2)" 2 where ap and eg are the semi-major axis and orbital eccentricity of star B, respectively, capable of generating spin-orbit misalignment through secular resonances 
primordially (aB ~ Õsa + wsb + wsc), and after the gas disk dissipates (weg ~ wsb + wsc, green shaded region). Typical PMS stellar rotation rates (1 day < Prot < 
10 days) are bracketed by dotted lines. The dashed line indicates the currently measured rotation period. The primordial misalignment shaded region is bracketed by 
Ma = 0.1 Mọ (upper line) and Ma = 0.005 Mọ (lower line) for rout = 15 au (red shaded region) and rout = 50 au (blue shaded region), where we require 
Tout > 3 Gp from tidal truncation of the disk by star B. The lower limit for the disk mass (Ma = 0.005 Mọ œ% 5 MJup) is consistent with the lowest gas-mass measurements 
of protostellar disks in the o Orionis Cluster (22). 

The post-formation shaded region is bracketed by R, = 1.511 Ro (upper line) and R, = 1.0 Rọ (lower line). We take 

My =5 Me: all other parameter values are taken from primordial and post-formation misalignment scenario sections. 
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Table S1. RV measurements taken during the transit of K2-290 Ac with HARPS-N@TNG 


Time (BJD) RV-19,760(ms—t) ørv (ms—?) 
2458599.511429 -15.5 9.0 
2458599.525665 -8.9 8.5 
2458599.539358 -7.4 8.1 
2458599.553687 -20.5 8.0 
2458599.665613 3.5 8.1 
2458599.680174 17.9 8.4 
2458599.694954 20.6 8.3 
2458599.708427 -3.7 8.5 
2458599.734446 16.1 8.1 
2458599.755443 -11.7 8.1 
2458653.461391 -25.3 5.6 
2458662.439687 -33.4 6.7 
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Table S2. RV measurements taken during the transit of K2-290 Ac with HDS@SUBARU. 


Time (BJD) RV+12,000 (ms!) ory (ms~?) 
2458646.904564 5.9 7.4 
2458646.915414 17.8 7.0 
2458646.926504 16.5 7A 

2458647.750561* -7.2 6.3 
2458647.761661* -34.5 6.3 
2458647.773090* 31.1 9.7 
2458647.784190 9.1 6.7 
2458647.795279 1.5 6.8 
2458647.808349 -13.4 6.5 
2458647.819449 -18.6 6.4 
2458647.830548 -18.0 6.4 
2458647.841638 -21.5 6.3 
2458647.852737 -12.8 6.2 
2458647.863827 -28.9 6.1 
2458647.874926 -18.8 5.9 
2458647.886016 -22.2 6.5 
2458647.897115 -25.7 6.6 
2458647.908215 -20.5 6.8 
2458647.921004 -14.5 6.7 
2458647.932104 -21:7 6.9 
2458647.943204 2.6 6.9 
2458647.954303 -11.2 6.9 
2458647.967083 2.2 6.7 
2458647.981662 -0.5 6.3 
2458647.996231 15.1 6.7 
2458648.009101 5.6 6.8 
2458648.020200 16.7 7.4 


Notes: *Excluded in the analysis due to them being affected by twilight or technical problems with the CCD read-out. 
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Table S3. RV measurements taken during the transit of K2-290 A b obtained from ESPRESSO@VLT. 


Time (BJD) RV-19,670(ms—!) ørv (ms—!) 
2458685.496193 -2.6 1.9 
2458685.507030 -2.4 1.9 
2458685.517732 2.4 2.0 
2458685.528853 -1.6 2.0 
2458685.539639 2.8 2.6 
2458685.550695 -0.1 2.6 
2458685.561700 -2.0 2.6 
2458685.572653 -5.1 2.1 
2458685.583376 -3.5 2.1 
2458685.594217 -1.8 2.1 
2458685.605165 -3.5 2.4 
2458685.621852 -3.9 1.8 
2458685.635662 0.1 1.9 
2458685.650111 23 23 

2458685.664976* -4.4 2.5 
2458685.679009* -6.1 2.2 
2458685.701177* -4.8 3.0 
2458685.728771* -5.7 2.1 
2458685.739713* -12.8 2.3 


Notes: *Excluded in the analysis due to the cloud coverage. 
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Table S4. K2-290 A system parameters obtained from transit data. 


We report median and 1-ø values. The “Prior” columns specify the distributions from which we draw our samples from; U 
denotes a uniform distribution, N (u, o) and T (u, c,a,b) denote a normal and truncated normal distribution, respectively, with 
mean p, standard deviation g, lower boundary a, and upper boundary b. 


Parameter Value Prior Value Prior 
Planet c Planet b 
Proj. spin-orbit angle A (deg) 153 £8 u 17375 u 
Proj. stellar velocity v sini, (km s7 +) 6.9408 N (u = 6.5,0 = 1.0) 6.9 + 0.5 N (u = 6.9, o = 0.6) 


Orbital period P (days) 48.3674 + 0.0003 u 9.2117 + 0.0002 u 

Time of mid-transit To (BJD) 2458019.1723 + 0.0003 u 2457994.77214 0: 001E u 

Scaled planetary radius Rp /R» 0.0684 + 0.0004 u 0.0196 + 0.0003 (i 

Scaled orbital distance a/R, 42.9 + 0.9 uU 13.1 0.6 u 

Orbital inclination i (deg) 89.34 + 0.06 u 88.2708 u 

RV semi-amplitude K, (m s71) 38.5177 N (u = 38.4,0 = 1.7) Lato 0.0, 0 = 2.2,a = 0,b = œ) 
Systemic velocity HARPS-N yna rPs-n (Mm s7!) 19, 76173 u E - 

Systemic velocity HDS yups (ms~+) —11,996+3 u - - 

Systemic velocity ESPRESSO yesrresso (m s7!) - - 19, 699.1493 

Linear limb darkening coefficient in Kepler-band u1, x 0.40 0.04 o =0.10,a = 0,b = 1) 0.44 + 0.08 T(u = 0.40, o = 0.04,a = 0,b = 1) 
Quadratic limb darkening coefficient in Kepler-band u2, x 0.104905 ao =0.10,a=0,b=1) O11 0-68 T(u = 0.10, o = 0.06, a = 0,b = 1) 
Linear limb darkening coefficient HARPS-N ui, HARPS—N 0.36 + 0.10 og =0.10,a=0,b=1) - - 

Quadratic limb darkening coefficient HARPS-N w2,HARPS—N 0.30 + 0.10 og =0.10,a=0,b=1) - - 

Linear limb darkening coefficient HDS ui,ups 0.444000 og =0.10,a=0,b=1) - - 

Quadratic limb darkening coefficient HDS w2,Hps 0.3640 og o =0.10,a =0,b = 1) - - 

Linear limb darkening coefficient ESPRESSO ui ESPRESSO - 0.37 + 0.10 T(p = 0.36, 0 = 0.10,a = 0,b = 1) 
Quadratic limb darkening coefficient ESPRESSO u2 nspresso - 0.30 + 0.10 T(p = 0.30, 0 = 0.10,a = 0,b = 1) 
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